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Graphene is a promising candidate for applications in spintronics. In this paper, Density Func- 
tional Theory method is used to calculate the band structure and magnetic properties of graphene 
on Ni(lll). Our results show that once there is anti-ferromagnetic order in graphene, an external 
electric field at the order of 10 9 V/m can induce a gap width difference of tens of meV for opposite 
spin states near the Fermi surface. 



I. INTRODUCTION 



Graphene is a two-dimensional material with many 
interesting propertied^! Behavior of Dirac Fermions 
in graphene has been extensively studied. Particularly, 
magnetic properties have attracted lots of interests. 

Graphene is considered a promising candidates for 
spintronics due to its long spin relaxation time and large 
spin relaxation lengtlP. Spin splitting is usually a de- 
sirable effect in potential applications. Although the in- 
trinsic spin orbit coupling (SOC) is quite wealP, con- 
siderable extrinsic SOC effects have been proposed and 
studiecP^, among which the Rashba SOC resulted from 
broken inversion symmetry received lots of attention. 
The Rashba Hamiltonian can be expressed as 



H = a(p x a) ■ E 



(1) 



a , p = fik , a and E respectively stands for Rashba 
strength parameter, electron momentum, Pauli spin op- 
erator and the electric filed. 

Graphene/Ni system is studied here for two reasons: 
One is that graphene films can be conveniently prepared 
on Ni substrate; the other is that Ni is a common can- 
didate for spin injection electrode and spin detection 
electrode. Published work suggests that the Rashba- 
induce d spin g ap can be as large as tens to hundreds 
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However, systematic study of the influences of Rashba 
effect on the Dirac fermions around the K point near the 
Fermi surface which greatly affect the transport proper- 
ties is still lacking. 

In this paper, we study few layer 
graphene(FLG)/Ni(lll) system using density-functional 
theory (DFT) by considering various factors including 
substrate interaction, layer number and SOC. 



II. METHODS 

For Firs t Principle Calculations LDA with PAW in 
VASP^EI] has been used to optimize the Few Layer 
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FIG. 1: (Color online) (a) Side view of bilayer 
graphene/Ni(lll) system. Three lower layers and two upper 
ones respectively stand for Ni and graphene. (b) Top view of 
top-fcc structure. 



Graphene (FLG)/Ni system, which according to previ- 
ous works is effective and sufficient^. The energy cut- 
off of optimization is 500eV and the convergence criteria 
is that the Hellmann-Feynman forces be less than 0.03 
eV/A for the geometry optimization. The graphene lay- 
ers and upper three layers of the total six layers Ni(lll) 
are allowed to relax. A vacuum slab more than 15A is 
applied to avoid interaction between supercells. 

The original lxl unit cell of Ni substrate is used, 
onto which the lxl unit cell of graphene is mapped. In 
monolayer case the graphene is of the well studied top- 
fcc site, as shown in FIG. 1(b); for bilayer the bottom 
lay er rem ains the same while the layer above is of hcp-fcc 
sit e! 20 * 21 !; in case of more than two layers of graphene, an 



2 



AB stacking model is used. This configuration is shown 
in FIG. 1(a). 

Then SOC is included into calculation to determine the 
influence of Rashba effect on the properties of monolayer 
case. 



III. RESULTS AND DISCUSSIONS 

The calculated band structures along K — T — K for 
1-4 layers of graphene on Ni(lll) are demonstrated in 
FIG. 2. We perform projected band structure calculation, 
in which the radius of red/blue spots demonstrates the 
weight of p z orbit of top layer graphene. FIG.2(a)/2(b) 
display band structure near Fermi surface for major- 
ity/minority spin of monolayer case, while that of bi- 
layer, trilayer and quadrilayer are respectively shown in 
FIG.2(c)/2(d), 2(e) /2(f) and 2(g)/2(h). 

As is shown in FIG. 2(a) and 2(b), the existence of 
Ni substrate seriously changes the band structure of 
graphene compared with freestanding situation with gap- 
less Dirac cone and linear energy dispersion. Due to 
the broken of the equivalence of the two sublattices of 
graphene caused by the substrate interaction, the Dirac 
Cone must be destroyed and a gap about 0.25eV opens 
up. The dense gray bands below the gap is those of Ni 
3d orbit. The ferromagnetism of Ni leads to a 0.4eV 
spin gap, indicating the spin-polarized characteristics of 
graphene layer. The detail magnetization is displayed 
in FIG. 3, in which red and blue circle represents mag- 
netism of C atoms parallel and anti-parallel to that of 
Ni atoms while the radius stands for the magnitude re- 
spectively. For convenience, the top layer Ni substrate is 
drawn as large gray circles. The graphene has an anti- 
ferromagnetic order, with one sublattice, namely A, un- 
der which is a first layer Ni atom, of magnetic momentum 
anti-parallel (O.QlfiB ) to that of Ni, and the other sub- 
lattice, namely B, of a parallel magnetism (0.03/iB ). 

For bilayer graphene, the parabolic energy dispersion 
is recovered since interation between graphene layers is 
weak. The Dirac point is lowered by about 0.4eV, sug- 
gesting strong charge transfer and vertical electric field, 
which can explain the small energy gap. And the spin 
gap is almost immeasurable, displaying the non spin- 
polarized character of top layer. Our calculation of the 
magnetization is also one magnitude lower than mono- 
layer case. In Fig. 2(e) to 2(f), the energy of Dirac cone 
is at 0.2eV. In quadrilayer case, the Dirac Cone is almost 
at the Fermi Level, as shown in FIG. 2(g) and 2(h). This 
can be understood since less charge will be transferred 
to upper layer. And we expectedly find no difference be- 
tween spins and negligible magnet moment. 

Thus from the discussion above, we conclude that the 
Ni substrate will result in a magnetization of graphene 
atoms and a vertical electric field. Both effects are 
weaken when layer number of graphene increase. Then 
Non-collinear Calculation is done for monolayer case 
since it has observable magnetization. We are inter- 
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FIG. 2: (Color online) Band structure near Fermi surface 
for majority /minority spin of monolayer (a) and (b), bilayer 
(c) and (d), trilayer (e) and (f), quadrilayer (g) and (h). The 
radius of red spots demonstrate the weight of p z orbit of top 
layer graphene. The insets show the detail near -K point and 
Fermi Level. 
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FIG. 3: (Color online) Calculated magnet moment of 
graphene atoms. Red/blue circle represents C atoms with 
magnetism parallel/anti-parallel to that of Ni atoms and the 
radius stands for the magnitude respectively. 



ested in the band structure along T — K near Fermi sur- 
face. The electric field E, as discussed above, is along 
the normal direction of Ni(lll) surface, the measured 
wavevector k is along r — to which the magnetism 
of the system is set perpendicular in plane. The calcu- 
lated band structure of monolayer case involving SOC 
is shown in FIG. 4(a), in which red/blue circles stand 
for bands of majority/minority spin with their radius 
representing weight of it orbit of graphene atoms. The 
kj — k symmetry of bands is destroyed due to the linear 
dependence of Rashba splitting on wavevector. To de- 
termine the effect on spin gap of Rashba and exchange 
splitting individually, we do calculation of the same sys- 
tem both with/without SOC. With SOC, the spin gap 
will be E — Eso + Eex, respectively representing the 
contribution of SOC and exchange splitting. If without 
SOC, only Eex is included. In this way, we can show 
the individual Rashba splitting in FIG. 4(b). One im- 
portant thing we should notice is that, for each spin in 
FIG. 4(a) there are two bands near the Fermi surface, one 
above the energy gap at K and the other below. Rashba 
splitting for upper/lower band is represented with grey 
squares/green circles in FIG. 4(b). It is critical that the 
Rashba splitting of upper (30meV) and lower (60meV) 
bands are of different sign, which means that the energy 
gap of one spin at K will be enlarged while the other 
shrinked. This can be understood since the upper and 
lower bands near K and -K are mainly contributed by 
sublattice A and B respectively, see FIG. 4(c) and 4(d). 
As has discussed above, the two sublattices are of oppo- 
site magnetism, which, according to Eq.(l), will result in 
opposite shift in spin gap. Due to the periodic boundary 
condition used in simulation, Eso decreases to zero at 
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FIG. 4: (Color online) (a) Calculated band structure of 
monolayer graphene /Ni(lll) system near Fermi level after in- 
troducing SOC. Red/blue circles stand for majority /minority 
spin bands, and their radius represents weight of graphene 
atoms' 7r orbits, (b) Rashba splitting of upper (grey square) 
and lower (green circle) band. (c)/(d) the contribution of n 
orbit of sublattice A/B to the bands, the radius indicates the 
weight. The insets show the detail near -K point and Fermi 
Level. 



the boundary of B^l Thus we evaluate it from near K 
and -K, which will be underestimated. The Rashba spin 
gap for upper band decrease sharply when getting away 
from K and -K(see the gray squares in FIG. 4(b)). This 
is because of the mixed contributions of both sublattices 
to the upper band at wavevector away from K and -K, as 
shown in FIG. 4(c) and 4(d). However, if we are focusing 
on band within 0.5eV around energy gap which mainly 
determains the transport property, the Rashba effect re- 
mains more than 30meV /-50meV. In cases of more than 
one layer of graphene, the Rashba effect on band struc- 
ture of top layer also becomes inconsiderable due to the 
almost disappearance of magnetism. 

In this case, we have shown the comparatively large 
Rashba induced extra spin gap for bands at K near Fermi 
surface. In addition, the anti-ferromagnetic order and 
two sublattice character of graphene lead to the discrep- 
ancy of energy gap between spins. One will be enlarged 
for about 45meV and the other shrinked for about the 
same magnitude. Here the electric field in Rashba Hamil- 
tonian is produced by the spontaneous charge transfer, 
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FIG. 5: (Color online) Schematic diagram to apply our 
method, excluding exchange splitting, (a) with zero elec- 
tric field in Rashba Hamiltonian, black lines stands for bands 
of both spin up and spin down, (b) with the electric field 
in Ni(lll) case, blue/red line respectively represents spin 
up/down, (c) with external electric field or strengthened mag- 
netization of graphene, the Rashba effect is enlarged and the 
system is converted from semiconductor to half metal. 

which is approximately 0.4eV/^P23. Such effect can be 
tuned through extrinsic electric field perpendicular to the 
interface, and an electric field strong enough may com- 
pletely close the gap for one spin and remarkably enlarge 
the gap for the other spin, which is the ideal situation 
in spintronics. FIG. 5 is a schematic diagram in which 
we exclude exchange splitting and only consider Rashba 
effect. FIG. 5(a) is the case with zero electric field and 
no Rashba effect, (b) is the Ni(lll) case. It may be 
converted from semiconductor to half metal if we futhcr 
apply external electric field or strengthen the magnetiza- 
tion of graphene, as shown in FIG. 5(c). Another critical 
implication is its application in valleytronics. Due to the 
linear dependence of Rashba splitting on wavevector, the 



degeneracy of the valley degree of freedom is broken and 
the shrink/enlarge effect is opposite at K/-K. This indi- 
cates that under proper configuration one spin may only 
transport at K and the other only at -K. Thus the Rashba 
effect we discuss here can also be an effective valley filter. 



IV. CONCLUSINOS 

In summary, graphene/Ni (111) are studied consider- 
ing several factors including magnetism induced by sub- 
strate and Rashba effect. At last, the Rashba effect is 
found to produce a 30 meV/-60meV extra spin gap on 
upper/lower band at K near Fermi surface, which can 
further adjusts the spin gap with Rashba effect. In addi- 
tion, the possibility to make discrepancy between energy 
gaps of different spin is discovered and interpreted, which 
provides a completely new method to efficiently optimize 
graphene/Ni system for the application in spintronics. 
And the opposite effect at K/K' makes it also a promising 
valley filter. This method is tunable and can be extended 
to other systems with similar characters. Actually, we 
expect such effect in all graphene systems with different 
magnetization of two sublattices. Graphene on semicon- 
ductor substrate is desirable for application and SiC is 
a promising candidate. Another perhaps better system 
is carbon nanostructures without any substrate, which 
can display antiferromagnetic behavior due to topologi- 
cal frustratiorP^. Graphene nanoribbon is an particularly 
expected choice, which is under way. 
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